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The Complex Scaling Method (CSM) provides scattering wave functions which regularize reso¬ 
nances and suggest a resolution of the identity in terms of such resonances, completed by the bound 
states and a smoothed continuum. But, in the case of inelastic scattering with many channels, the 
existence of such a resolution under complex scaling is still debated. Taking advantage of results 
obtained earlier for the two channel case, this paper proposes a representation in which the conver¬ 
gence of a resolution of the identity can be more easily tested. The representation is valid for any 
finite number of coupled channels for inelastic scattering without rearrangement. 


I. INTRODUCTION, NOTATIONS 


As is well known, the CSM converts the description of resonances by non-integrable Gamow states into one by 
square integrable states while leaving the discrete spectrum unchanged [1], Cuts describing the continuum are rotated, 
however, but this may be advantageous, since they are thus disentangled when their thresholds differ from one another. 
(We are not interested, in this paper, in the case of channels with identical thresholds.) It is then expected that the 
continuum corresponding to such rotated cuts makes a much smoother contribution to the calculation of collision 
amplitudes, level densities, strength functions and sum rules [2] [3], since narrow resonant processes have been assumed 
to be peeled out explicitly by the CSM. The CSM Hamiltonian, unfortunately, is not hermitian any more, and it is not 
obvious that a resolution of the identity in terms of its bound states, resonances and presumably damped continuum 
is possible. For the one channel case, convincing arguments have been advanced a long time ago [4] to prove that this 
resolution exists. More recently [5], a detailed investigation of the case of two channels, coupled by straightforward 
potentials, generated a contour integration of the usual Green’s function which provided the identity resolution. The 
task was made reasonably easy by the small complication of the Riemann surface in that case. The purpose of the 
present paper is to capitalize on the methods used for that two channel case and attempt a generalization to any 
finite number of channels, despite the more complicated nature of the relevant Riemann surface. We shall assume, 
naturally, that there already exists, derived from single poles and usual cuts, a resolution of the identity for the initial 
Hamiltonian, before its modification by complex scaling. Our problematics would be meaningless otherwise. 

Several earlier studies, in particular by [6] [7], have been concerned with a description of resonances with square 
integrable states, without complex scaling. They did not restrict to the consideration of just simple poles of the 
S'-matrix and investigated how one might, as rigorously as possible, define initial wave packets for the description of 
decaying states; the non purely exponential nature of their decays received a detailed attention, via the analysis of 
their time dependent evolutions. The present paper, however, will be content with a Gamow definition of resonances, 
by means of simple poles; our aim is just to generate a resolution of the identity, with time independent states 
extending to asymptotic regions. For earlier searches of a complete basis of states, including resonances, but within 
a compact interaction volume, we may refer to the review by [8] of /?-matrix methods and their extensions. See also 
[9] and in particular the comparison of “class B” and “class D” theories. 

In this paper, we shall again assume that all potentials Vj n (r) driving the channels and their couplings are local 
and so short ranged, Gaussian-like for instance, that the 2 N Jost solutions of the N coupled equation system, 
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exist and are analytical in the whole complex domain of all the momenta kj. The radius r runs from 0 to +oo, 
obviously, and the number N of channels is taken as finite. As an additional technicality we also assume, naturally, 
that the products V ln ip n j do not diverge for r —► 0 when singular solutions of Eqs.(l) are considered. 

We select the threshold of the lowest channel as the origin of the complex energy plane, hence E = k 2 . The 
other channels with their physical thresholds E*, which are real and positive numbers, now define channel momenta 
according to, Ej = k 2 = E — e he E*. Notice that, given a real number E* defining a physical threshold, the usual 
complex scaling where p 2 becomes e~ 2lB p 2 and r becomes e l0 r does not change E* and rotates the corresponding cut 
by angle —26. But here, we have a slightly different representation, because the Hamiltonian has been multiplied by 
e 2lB . Hence kinetic operators in our Hamiltonian H , see Eqs.(l), are just —d 2 /dr 2 , every cut rotates back into being 
“horizontal” and starts from e 2l0 E*. For time dependent studies, it will make sense to scale time, conjugate of energy, 
by a factor e~ 2lS . This will prevent those resonant wave packets, the energies of which have a positive imaginary part 
as eigenvalues of H, from exploding when t —> +oo. 

Also in this paper no rearrangement is allowed, channels are defined by just internal excitations of the projectile 
and/or the target, hence all reduced masses are equal. Finally we exclude from this paper the consideration of 
abnormal thresholds; we shall only discuss the case of “square root thresholds”. This is generic enough. 

It is understood here and from now on that a first subscript, such as i or n, denotes the component of each wave ip 
in channel i or n, then that any superscript, ±, or second subscript, j, denotes the boundary condition which defines 
ip. For a Jost solution /// the boundary condition that we choose is “asymptotic flux eih^-tW 2 ) in channel j and 
no asymptotic flux in the other channels”. It is well known that for r —> 0, the components of such Jost solutions are 
proportional to (kir)~ li ( 2f; — 1)!!. For a regular solution cpj, the boundary condition that we choose sits at r = 0 and 
reads, “lim r ^ 0 Pij(r) = 0 V* ^ j, while, for i = j, then lim r ^ 0 <J>jj{ r ) = 1/(2 lj + 1)!!. 

Following Newton [10], it is convenient, given E and r, to set the column vectors ip.j into a matrix <f?(E,r) of 
regular solutions and the Jost solutions /'t (resp. /“) into a similar matrix i + (E,r) (resp. f _ ). It is also convenient 
to notice that <1?, viewed as a function of the kj' s as if these were independent momenta, is even under any reversal of 
a kj into —kj. Such is not the case for f + ; analytic continuations in either energy or momenta planes can introduce 
one (or several) /“’s into f + . 

For our oncoming argument we must use the Wronskian matrix with matrix elements the Wronskians W (/+,, p. n ) 
of the Jost solutions /+, with the regular ones p. n . This, for s waves, is the transposed of f+ at r = 0, 

W(E) = ~f+(E,0), (2) 

and for other angular momenta is only a slight modification of f + (F,0). (Rather than just f + (F,0) one must use 
limits of products (fc,;r)^ i /l'/(2f/ — 1)!! at r = 0, explicitly, but we will disregard this technicality.) The Green’s 
function G is then found as, 

G(E, r, r') = &(E, r) [W(F )] _1 i + {E, r') if r < r', G(E,r,r') = f + (E,r) [W (£ , )] _1 &(E,r') if r >/. (3) 

Here each tilde ~ means transposition; we refer to [10] or to Appendix A of [5] for the derivation of such formulae 
for G. Despite different formulae whether r > r' or r < r ', and the lack of hermiticity, G is symmetric, namely 
G(r,r') = G(r',r). 

It will be noticed that the CSM, as we describe it by the system of Eqs.(l), locates thresholds on a segment of the 
complex E plane with slope 26, extending from E = 0 to e 2lS E )/, and that the channel cuts are rotated back into 
being “horizontal”. Conversely, bound states lie on a negative semiaxis rotated by 26 and resonances are rotated by 
26 as well. This slight change of representation changes nothing to the physics, obviously. For trivial technical reasons 
[5], we normalize energy units so that E^ = 4. Also we shall use a short notation, k = ki and K = kw- We show in 
Figure 1 the cut energy plane in an illustrative, four channel situation when 6 = 7r/6, E% = 1.5 and E£ = 3.5. 

Equipped with this slightly unwieldy formalism, we can now investigate whether there exists a representation, and 
an integration contour, such that the traditional integral, T = J dE G(E, r, r 1 ), calculated in two different ways, 
generates a resolution of the identity. This question of a representation and a contour is the subject of Section II, the 
main part of our argument. Additional considerations on the two ways of calculating this integral make the subject 
of Section III. A discussion and conclusion are proposed in Section IV. 




FIG. 1. P-plane. Physical cuts for a four channel case when 0 = n/6, E% = 1.5, E% = 3.5 and E% = 4. Lowest channel, 
heavy full lines, highest channel, heavy dashed lines, intermediate channels, lighter full lines. The dotted segment with slope 
7t/3 is the locus of thresholds (big dots) in this representation. 


II. REPRESENTATIONS AND CONTOURS 
A. Energy plane 


From Fig. 1 it is intuitive that one could start, for instance, from +oo along the lower rim of the lowest channel 
cut, return to the origin, E = 0, proceed to +oo again on the upper rim, then join there the lower rim of the second 
cut, return to the threshold of this second cut, go to e 2lS E% + oo along the upper rim, join the third cut lower rim at 
infinity, etc., until arriving at e 2lS + oo along the upper rim of the highest channel. Then the contour would be 
closing at infinity by means of an almost complete circle, counterclockwise, terminating at the starting point, namely 
at +oo on the lower rim of the lowest channel. 

Along such a contour, it would be necessary to investigate the behaviors of the ingredients f + , W and $ of G. 
Furthermore, information is needed about the singularities of G inside the contour; indeed, residues of simple poles 
are essential for a calculation of f dEG(E) by Cauchy’s theorem; one also needs reasons why no singularities higher 
than simple poles occur. 

The representation discussed in the next subsection makes easier the needed investigation, for it opens two of the 
cuts and limits the discussion to situations where all momenta have semipositive imaginary parts, Qkj > 0. 


B. Pseudomomentum plane 


A generalization from [5], where there were two channels only, the present “P representation” consists in joining 
the upper rim of the lowest cut and the lower rim of the highest cut, and in opening both cuts, by rational formulae, 
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k = P + Q 2 /P, K = P-Q 2 /P, 

where Q = e l6 makes a short notation for our scaling of energies such that E^ = 4 and k 2 — K 2 = 4 Q 2 . Trivially, 
P is the average ( k + K )/2 of k and K. The point is, despite an obvious failure to open additional cuts, P also give 
the “dominant” part of any other momentum when 3P —> +oo. Indeed, when |P| is large, say |P| >> 2, then an 
asymptotic value can be defined for kj, j ^ 1 , j ^ N , according to the rule, 

kj = ( k 2 - Q 2 E = (P 2 + 2 Q 2 - Q 2 E* + Q*/P 2 )? =P + Q 2 ( 1 - E*/2)/P + 0{P~ 2 ). (5) 

Thus the semicircle at infinity in the upper P plane corresponds to Qrfcj > 0, Vj. This is of critical value for the zoology 
of our Jost functions and it is expected that this semicircle properly closes the integration contour under design. 

Set now P = x + iy and short notations c = cos 29 and s = sin 29. A trivial calculation separates the real and 
imaginary parts of the (complex) energies driving each channel, 

(x 2 +y 2 ) 2 $l(k 2 ) = [(a ; 2 + y 2 + s)(x + y) + (x - y)c] [(a : 2 + y 2 - s)(x - y) + {x + y)c) - E*(x 2 +y 2 ) 2 c, (6) 

and 


(a : 2 + y 2 ) 2 3(fc 2 ) = 2[(a : 2 + y 2 )x + xc+ys] [(x 2 + y 2 )y + a :s - yc ] - P](a ; 2 + y 2 ) 2 s. (7) 

and it is trivial to recover the images, in this new representation, of the cuts displayed in Fig. 1. Polar coordinates, 
with P = pe lrl , can be also be used to decribe the j-th cut from Eq.(7) by, 
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p 2 sin 2 77 H-^-= (Ea — 2) sin 29. 
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FIG. 2. P plane. Cuts for the same four channel case, 9 = 7 r/ 6 , = 1.5, P 3 = 3.5 and EX = 4. Opened cut for lowest 

channel, heavy full line. Opened cut for highest channel, heavy dashed line. Intermediate channel cuts, not open, lighter full 
lines. The dotted segment is the locus of thresholds (big dots) in this P representation. 



Results are shown in Figure 2 for the same special case as Fig. 1. As in [5], the lowest channel is represented by 
the heavy, shoulder shaped line, that starts from —oo on the real P axis, bends up, then backs into the origin P = 0, 
where it terminates with a slope 29. Along the curve, k is real and runs from — oo to +oo, covering both rims of the 
initial cut. The threshold k = 0 is represented by P = iQ = e*^ + 2 7r ). Partner points where k —k obtain under the 
symmetric transformation P —Q 2 /P. In the same way, for the highest channel, K runs with real values along the 
heavy dashed line, from — oo at P = 0 to +oo at the end of the positive 3?P semiaxis, via K = 0 for P = Q. The 
transform, P *-> Q 2 /P, makes partners with opposite values of K. 

The other cuts remain cuts. Their thresholds lie on the image, shown as a dotted line again, of the segment already 
pointed out at the stage of Fig. 1. Because both <! ft(k 2 ) and c A(k 2 ) vanish for such points, it is easy to eliminate E* 
between the right hand sides of Eqs.(6,7) and obtain the condition for such a locus, 

x 2 + y 2 = 1, (9) 

a very simple result indeed. With |P| = 1, the positions of the thresholds are easy to obtain. The special cases j = 1 
and j = N give the argument 77 = ArgP as 77 = 9 + 7 r /2 and 77 = 9 , respectively. This was already known from [5]. 
The function sin 2?7 + sin(40 — 2if), see Eq.( 8 ), decreases monotonically when 77 increases from 9 to 9 + 7 r/ 2 , hence a 
unique solution for each E*, and an obvious symmetry about 9 + 7 r /4 corresponding to the symmetry about E* = 2. 
Then each intermediate cut generates, from Eq.(7), an image which joins its threshold to the origin P = 0, while kj, 
a real number along this image, runs from 0 to ± 00 , according to the rim. The image lies between the heavy full and 
dashed lines, and, being pinched between them at P = 0, also reaches the origin with slope 29. While the pinching 
makes numerics slightly difficult, it is easy to verify analytically from Eqs.(6,7) that infinitesimally away from both 
rims of such an intermediate cut, but inside the wedge created by the heavy line curves, ‘Akj remains positive. 
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FIG. 3. P plane. Again 9 = n/6. Cut for the channel defined by E 2 = 1.5. The center line, between dots, is the cut. Cut 
then continued for negative energies in the channel. Additional lines, lower rim (leftmost curve) and upper rim (rightmost 
curve), respectively. Both rims extended below threshold. Heavy line bar, connection between extended rims. 

To illustrate our full control of the various 3/cj’s provided by this P representation, whether inside the wedge or near 
the positive infinity semicircle, we show in Figure 3 the cut corresponding to E%, and its continuation beyond threshold. 
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“semicircle”. Simultaneously, we generate rims of the cut, and beyond again below threshold. To generate rims, we 
use Eq.(7), or as well Eq.(8), with E% replaced by E | — 0.2 and E% + 0.2 for the lower and upper rim, respectively. 
(The choice ±0.2 was made for graphical convenience, but we tested much smaller intervals, naturally.) The dots 
represent P = 0, where the channel energy is infinite, and the threshold, where it vanishes by definition. Like the cut, 
the rims are pinched by the wedge. 

Then we show in Figure 4 the trajectory of when P follows this cut from P = 0, to the threshold and beyond. 
Notice that, E 2 being real along the line, then the imaginary part of £3 = E 2 + e 2ie (E% — £ 3 ) is obviously negative. 
This does not prevent a choice of £3 with 9^3 > 0, generating the leftmost trajectory in Fig. 4. Simultaneously, we 
show the trajectories of k 2 from both rims of the same cut. The left hand side (when seen in Fig. 3) rim induces 
$tk 2 —> —oo when P —> 0, with an infinitesimally positive 3 k 2 . Conversely the right hand side rim induces Wt 2 —> +oo 
when P —> 0, with still an infinitesimally positive $ sk 2 . When we go from either rim towards the upper semicircle at 
infinity, this induces Sfc 2 —> +oo, as expected. The rims can be connected by any small path, see the bar above the 
threshold in Fig. 3, and the values of k 2 along the rims can be smoothly matched, see the curved bar in Fig. 4, the 
trajectory of k 2 when P follows the bar in Fig. 3. Generalizations to every kj in every part of the wedge are trivial. 
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FIG. 4. k 2 , kz planes. Still 9 = 7 r/ 6 , E 2 = 1.5 and £3 = 3.5. Leftmost curve, trajectory of k$ when P follows the central line 
of Fig. 3. Intermediate curve, trajectory of k 2 for extended lower rim, see leftmost curve in Fig. 3. Rightmost curve, trajectory 
of k 2 induced by extended upper rim, see rightmost curve in Fig. 3. Heavy line curved bar, connection trajectory for k 2 when 
P turns around the threshold, below it. 


C. Contour 


To synthetize this Section, the P representation defines a physical sheet similar to the physical sheet of the energy 
plane. The region of interest is that region above the two curves which open the cuts for the lowest and the highest 
channels, while cuts remain for the intermediate channels. All momenta inside the wedge, and all the way to the 
upper semicircle at infinity, can be defined with positive imaginary parts. A contour can be found, following all cuts 
and closing at infinity in the upper plane. 

The intuition which was present in the E representation can be substantiated in the P plane. Start from —oo on the 



real axis, follow the “opener curve” which corresponds to the lowest channel, all the way to P = 0. From there, follow 
the lower rim of the cut corresponding to the second channel, back to its threshold, then turn around the threshold to 
follow its upper rim, down to P = 0. In turn, follow the lower rim of each intermediate channel, then its upper rim. 
After bouncing N — 1 times at P = 0, follow the “opener curve” corresponding to the upper channel, until P —> +oo 
on the real axis. Then close the contour by means of the upper semicircle at infinity. In the next Section, we shall 
investigate what happens to the integral, T = J dE G (E, r, r'), when considered along this contour in the P plane. 


III. THREE CONTRIBUTIONS TO THE GREEN’S FUNCTION INTEGRAL 


A. Upper semicircle 


At infinity in this upper P plane, the integration weight, dE = 2 (P — Q 4 /P 3 ) dP, boils down to 2 P dP. All the 
TV distinct Jost solutions boil down to exp [?'(Pr — |^ 7 r)] in their respective “flux channel j”, while vanishing in the 
other channels. At the same time, the TV distinct regular solutions similarly boil down to sin (Pr — n) /P in their 
respective flux channel and vanish in the other channels. The Wronskian matrix boils down to the TV-dimensional 
unit matrix. 

Assume now r > r', for instance, and thus consider the second of Eqs.(3). The product f + [W ]^ 1 3> boils down to 
a diagonal matrix. Its j-th diagonal element reads, 

[ 2 dP e i ( Pr- ^' 7r / 2 ) S i n f Pr ’ _ . ( 10 ) 

and can be easily calculated by reducing the semicircle back to the real P axis. The result does not depend on j, 
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dP e i(Pr -3^) 


exp (iPr — TF,—^ ~ ex P (*^ 7 ^ — iPr' 


= 2 ?' 7 r[( 5 (r + r') — S(r — r') 


( 11 ) 


It is trivial to verify that the same result is obtained if r < r'. Furthermore the term 5(r + r') cancels out in the space 
of regular radial waves. Hence the contribution X sc of the semicircle makes nothing but the multichannel identity, 
multiplied by (— 2 * 7 r). Notice that, differing from [2], this identity is not multiplied by a factor depending on 9 , since 
for us the ends of the semicircle, —00 and + 00 , both lie on the real P axis. 


B. Continuum 

It makes no difference here whether we consider the contribution of one of the “opener line” or that of one of the 
intermediate cuts. For in both cases we group partner terms. Such partners either come from a transform P ±Q 2 /P 
or from opposite rims of the intermediate cut under consideration. What is important to notice is that momenta 
retain their finite and positive imaginary parts and do not change when we compare two partner points, except that 
momentum specific to the opener line or the cut. For that momentum, which is real, “partnership” means kj <-> —kj, 
with still an infinitesimal positive imaginary part. Keeping in mind that $ is even under such a momentum flip, the 
contribution of such a continuum thus reads, if r > r' for instance, 

nOO 

Xj = / 2 kj dkj D j(E, r ) 4>(P, r'), (12) 

Jo 

where Dj (P, r) represents the following difference between partners, 

D j(E, r ) = f +(E,r) [W(P )]” 1 - f + (-kj,r) [W (-kj)]- 1 , (13) 

a discontinuity across the cut. The notation used here takes advantage of the fact that dE = 2kj dkj, and that kj is a 
convenient label along the line or the cut. The first term, f + (P, r) [W(P)] -1 , in the right hand side of Eq.(13) clearly 
comes from the upper rim. The notation that we use for the second term, f + (— kj,r) [W(— kj)] -1 , indicates that, 
because of analytic continuation in the physical sheet around the threshold, one Jost solution f~ now makes the j’-th 
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be a rank one dyadic. An elementary proof of this dyadic result was given in Appendix C of [5]. Nothing changes in 
the argument if r < r'. 

As a consequence of the dyadic nature of D^, and of the symmetry G(E,r,r') = G(P,r',r), hence of the same 
symmetry for discontinuities across cuts, there exists as a column vector a solution <f> M j of Eqs.(l) that is able to 
represent symmetrically both D j(E,r) &(E, r') and <J?(P,r) D j(E,r') in a self dual way as an outer product, 


*i = 


2kj dkj 


<l>.j(E,r) 4>.j(E, r') 
V{E) 


(14) 


This solution belongs to the set of regular solutions, naturally, because of the regularity of G at both r = 0, and 
r' = 0, illustrated by the presence of $ in Eqs.(3). The exact natures of this <j)j and of the “normalizing” denominator 
V are discussed in the Appendix. 


At this stage, the full integral along the full contour thus gives the sum of the multichannel identity and “pseudo¬ 
projectors on the continuum”, one pseudoprojector for each channel, 
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dE G(E, r, r’) 
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^.j(E,r) <j.j(E,r') 
V{E) 


(15) 


The next subsection shows what happens if the same integral is evaluated by means of the Cauchy theorem. 


C. Residues at poles 


We assumed that, before complex scaling, namely for 9 = 0, there existed an identity resolution in terms of unsealed 
bound states and unsealed scattering states. In other words we assumed that the corresponding, unsealed G (P) shows 
only isolated, simple poles, besides the physical cuts. Such poles can be on the real E axis of the physical sheet, 
describing bound states, or away from this axis, then describing resonances or antiresonances. The point is, now, that 
the CSM cannot change the nature of such poles [1]. Within our description by Eqs.(l), the CSM just rotates such 
poles by 29 in the energy representation, along circular arcs, concentric around E = 0. In the P representation, the 
images of such arcs are also concentric arcs, with angular extension 9 only. This is trivially seen from the equation 
which, for each initial position £ of a pole, defines those values of P which represent e 2l9 s, 

(P + e 2ie /P) 2 = e 2i0 e. (16) 

Indeed, 9 disappears from this equation if one sets P = e* e Po, where Pq solves for the initial position e. It can be 
concluded that only simple poles will be found when a finite 9 is used for our CSM. Notice, incidentally, that for e 
real and negative (bound states), the P representation will align poles along the axis with polar angle 9 + 7r/2, further 
than the circle with radius 1 that we found as the locus of thresholds. There will be no such alignment for resonances. 

For the calculation of T by Cauchy’s theorem, poles are not due to either f or <I>, since these, as functions of E 
or P, are regular. Only the divergence of W^ 1 can create poles. The situations of interest are those when the roots 
of the determinant, det W, are located inside the integration contour. We know that such is the case for the bound 
states. Depending upon 9, some resonances may also rotate into the domain. It is already known that only simple, 
isolated poles occur. The only question to solve is, what is the residue of G at such a pole. 

Residues of G at its poles will now be obtained from derivatives d/dE. That is equivalent to a calculation in the P 
representation, anyhow, and slighly easier. We shall use short notations in which the dependence of 3?, f + , W, upon 
r, and/or r' and/or E will be most often understood. However, at those energies E v where a pole occurs, we use an 
explicit subscript v to specify that such quantities <F, ... , W are evaluated at E v . 

Poles occur because of W -1 . Hence, we must only find the residue, 

Kv= lim {E-E^W-^E), (17) 

E — >E l/ 
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At a (simple!) root E„ of det W(E), there is necessarily one, and just one, null right eigenvector A„ of W. Similarly 
there is one, and just one, null left eigenvector A' v . We write them as columns and normalize them by the condition, 

K A„ = I- (18) 

Then the divergent part of W -1 in a neighborhood of E v is nothing but the truncation, 


W 


-l 

tr 


A'„W (E) A v 


(19) 


where there is an explicit dependence on E in the denominator. This denominator, a number, vanishes at E = E v . 
As a matrix element of W it is nothing but the Wronskian of the following two waves, F = f + A(, and £ = <E> A„. The 
former, F, is irregular, the latter, £, is regular. While A„ and A' v do not depend on E, since they were defined at 
E = E v , both F and £ depend on E, via f + and <l>. When their Wronskian vanishes, F and £ become proportional to 
each other, and there exits a number c such that F v = c£„. This special wave is both a mixture of regular solutions 
and a mixture of Jost solutions, with positive imaginary parts in the momenta driving all Jost solutions. Therefore 
it decreases exponentially in all channels when r —> oo and it is square integrable as well as regular. As expected it 
represents either a bound state or a regularized resonance. 

According to Eqs.(17,19), the residue under study comes from just the reciprocal of the derivative of the Wronskian 
of F and £, 


1Z U = 


KK 


KW(E) A„ /dE \ E = Ev 


( 20 ) 


In short, we must calculate the derivative of a Wronskian with respect to the energy, d |A(, W(E) A^J /dE. To help 
manipulations with Wronskians, define an operator matrix U with matrix elements the CSM potentials, completed 
by the centrifugal barriers and the thresholds, 


Uij - e 2W Uij (,'"/■) + 6.. 


e 2if> E* 


tj{tj + 1 ) 


( 21 ) 


Then elementary, but slightly tedious manipulations, which are already described in [10] or in Appendix B of [5], give 
the remarquably simple result, 


A/W(E)A V /dE \ E=Ev = - c / dr£(E v ,r)Z(E v ,r). 


( 22 ) 


Then the constant c cancels out between this and the numerators of f+ 1Z V and 1Z V f /~, which make the same, 
symmetric formula anyway, whether r > r' or r < r', since F v = d/ v . 

Summing upon all such residues obtained at roots E v of det W above the “opener” curves in the P upper half-plane, 
the contour integral reads, 


I(r, r') 


V 


*(E v ,r) A„ A„ *(E v y) 

/ 0 °° dr" A v 4>(E u ,r") *{E v ,r") A, ' 


(23) 


Here we state again that the column vector A„ is the null, right-hand side eigenvector of W (E„), namely W (E„) A v = 
0, then the column vector <&(£),) A^ of wave functions is the wave function of the bound state or resonance, and the 
denominator plays the role of a “Euclidean-like square norm”. This denominator is non vanishing; this corresponds 
to the hypothesis of single, isolated poles. All these are labeled by v, a discrete index, or as well by P v , an isolated 
root of W if viewed as a function of P. 


D. Completeness 


Since the three contributions I s . 




and I are obviously related by T sc + ]>A 2j- = T. it is trivial to equate sc , 
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and y^. , Ij, the latter term making the pseudoprojector upon the continuum for all channels. Naturally, in practical 
calculations, a cutoff and some amount of discretization will be necessary to integrate such continuum terms, but the 
P representation provides a suitable frame for testing the convergence of such a resolution for sum rules, level densities 
and similar observables. Notice that, because of the use of complex, self dual bras and kets in the resolution, such 
cutoff and discretization manipulations may generate spurious imaginary parts for the expectation values of hermitian 
observables. For a discussion and possible interpretation of imaginary parts in individual matrix elements, we refer 
to [11]. But, when summed upon all discrete and integral terms provided by the resolution, such imaginary parts 
must add up to a negligible, spurious noise compared to the real parts. This requested cancellation makes one more 
criterion to validate numerical operations. 


IV. DISCUSSION AND CONCLUSION 


Once again we used the ABC theorems [1] to locate the discrete spectrum at trivially rotated positions deduced 
from the discrete spectrum of an initial, hermitian Hamiltonian. The topological similitude provided by the CSM 
rotation warrants that, as long as there are no double poles or higher singularities with the initial Hamiltonian, the 
same will be true with the CSM Hamiltonian. 

Then it was not very difficult to find a representation which allows a suitable contour integration of the Green’s 
function. There was still a slightly complicated Riemann surface to handle, for the number of cuts was reduced to 
N — 2 only [12], but we took great care, including a few numerical, illustrative examples, to show that all cuts in the 
new representation are well understood, all thresholds are easily located, all complex momenta to be used for proofs 
have positive imaginary parts in a physical domain of a suitable sheet, and in general that all technicalities are sound. 

This proof of the CSM completeness for N channels is restricted to a finite number of well separated channels, 
normal square root threshold singularities, in a purely inelastic situation, without rearrangement, and with short 
ranged forces. The case of long range forces makes a more difficult question, indeed [13] [14]. But our restrictions 
still allow a large class of practical problems, and for instance in nuclear physics, a very large number of collective 
resonances can be described by the coupled channel equations that we studied. 
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Appendix 


We give here in some detail a description of that regular solution (j)_j which accounts for the discontinuity of the 
Green’s function across a cut. For the sake of pedagogy, we set the channel number to be N = 4 and shall consider 
only what happens for, e.g., the second cut. Generalizations are obvious and left as an exercise for the interested 
reader. In a condensed notation, we write the upper rim Wronskian matrix as, 
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(24) 


where, for instance, b is the Wronskian of f~[ with ip ,2 and o is the Wronskian of f \ with ^. 3 . The inverse of W 
reads, trivially, 


(25) 


where det u is the determinant of W„ and the prime symbols denote the corresponding cofactors. For the lower rim 
of the second cut a substitution occurs for the second row of W„, hence the lower rim Wronskian matrix reads, 


W, = 


a b c d~\ 
q r s t 
i j k l 

m n o p J 


(26) 


where, for instance, t is the Wronskian of / 2 with <£. 4 . Accordingly the inverse matrix becomes, 


Wf 1 = {deUY 


r j" 

' k" 
Id” h! I” 


c" ff 


p" J 


(27) 


where doubleprime symbols denote new cofactors, but the cofactors of { q , r, s, f} are the same as those of {e, /, g, h}. 
Again with a transparent, condensed notation, we set, for the upper and lower rim, respectively, 


(28) 


f 22:/32;/42}> while {Q.R,S,T} = 
{/ 12 , f '22 , / 32 > F°r r < r ', the discontinuity to be studied corresponds to the transposed of Eq.(13), and reads, in 
a condensed notation, 

D 2 (r') = Wf 1 P(r') - Wf 1 f ; (r'), (29) 

The subscript 2 for the cut and the r' dependence will be now understood and we shall use trivial identities to analyze 
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D = [Wf 1 + Wf 1 (AW) Wf 1 ] f“ - Wf 1 (f“ + Af) = Wf 1 (AW) Wf 1 f u - Wf 1 Af, 


(30) 


where AW = W/ — W u and Af = f ( — f . The point is, both modifications A are just substitutions for second rows; 
they boil down to dyadics, 

<g> [<7 — e r — f s — g t — h], Af = 


AW = 


® [Q-E R-F S-G T-H]. (31) 


(Our use of the tensor product symbol 0 is actually superfluous; we just want to stress the matrix product of a column 
by a row.) The next point is, then, that a global dyadic form for D emerges, 
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Furthermore, from the very definition of matrix inversion, we see that 
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hence 
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and D simplifies into 
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For r <r' the complete discontinuity 4>(r) D(r') of G(r, r') thus reads 

deti det u 3>(r) D(r') = <j>(r)E(r'), 


with 


4 > = e' ip. i + f ip. 2 + g' ip. 3 + h' ip a, 


(33) 
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and 


S = (qa + rb' + sc + td') f \ + {qe + rf + sg' + th') f% + ... + (, qm' + rn + so' + tp) - det u f 2 . (38) 

Both (f> and S are column vectors and relate to the second cut, hence they should actually read <f >.2 and 3.2 in a 
notation compatible with Eqs.(14,15). We omitted such subscripts, for the sake of conciseness. 

It may be convenient to take advantage of the cofactor nature of all the coefficients a ', ... p'. This gives indeed the 
formal, but condensed formula, 
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(40) 


(41) 


can themselves be interpreted, after keeping track of signs, as cofactors for the last column of the determinant shown 



In the 27V-dimensional space of solutions, it is known that the 2TV Jost solutions and the N regular ones are related 
by a formula such as, 


f = $ W ' 1 w + f+ w ~ 1 W W ' 1 w, 


(42) 


where W is the same as W u , while W_ is the analog of W if one replaces each /+, by its partner f7 rn . Then w is a 
diagonal matrix, defined from the Wronskians = —2ik m 5 mn . It will be noticed from Eq.(42) that, if we 

expand an f~ t on the basis spanned by all the p. n and all the /+ , the regular components of f~ n are provided by 
the TO-th column of the matrix product W ^ 1 w. 


It is known that 5 always belongs to the subspace of N regular solutions. In our illustrative example where N = 4 
and we studied the second cut, our H, according to Eq.(38), is a superposition of five solutions, namely all the /+ and 
one f~ only, /if. After an expansion of f7 2 on the basis spanned by the ip, n and the /+, all its irregular components 
must cancel out those preexisting irregular components of E seen from Eq.(38). (For the sake of rigor, we verified, by 
brute force calculations when N = 2, 3 and 4, that the components /+ do vanish out.) Thus we may consider the 
regular components only, coming from just 

The weight of fj 2 is, according to Eq.(38), — det u ■ We must therefore find the second column of, 
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(43) 


(44) 


The generalization, S.j = 2 ikj 4>,j is obvious. For any channel number N and any j-th cut, both (f) and 5 correspond 
to the j -th column of W , hence to the cofactors of the j-th row of W. There is no need here to specify W„ or W;, 
because the relevant cofactors are the same on both rims of the cut. The fact that H and (j> are the same except for 
the factor 2 ikj gives the same, symmetrical result whether r is larger or smaller than r'. And the denominator present 
in Eqs.(14,15) reads, when all factors are collected, 
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